Magnetic structure and phase diagram in a spin-chain system: Ca 3 Co 2 06 
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The low-temperature structure of the frustrated spin-chain compound Ca3Co2C>6 is determined 
by the ground state of the 2D Ising model on the triangular lattice. At high-temperatures it 
transforms to the honeycomb magnetic structure. It is shown that the crossover between the two 
magnetic structures at 12 K arises from the entropy accumulated in the disordered chains. This 
interpretation is in an agreement with the experimental data. General rules for for the phase diagram 
of frustrated Ising chain compounds are formulated. 

PACS numbers: 75.25.+Z, 75.30.Kz, 75.50.Ee 



A combination of a low dimensionality and frustration 
in magnetic compounds, which contain weakly bound 
Ising spin-chains packed into a two-dimensional (2D) 
frustrated lattice, gives rise to a complex magnetic be- 
havior. Antiferromagnetic (AFM) Ising chains form the 
triangular lattice in well-known Ising-chain compounds 

CsCoCl 3 , CsCoBr 3 0, 0, 0, S]. In these substances 
a partial AFM interchain order (honeycomb magnetic 
structure) appears just below the ordering tempera- 
ture. While the temperature decreases, two magnetic 
phase transitions occur to low-temperature ferrimagnetic 
phases. A theoretical explanation of the magnetic struc- 
ture of CSC0CI3, CsCoBr 3 was proposed by Mekata in 
the framework of a mean-field theory [l[. To describe 
the low-temperature phases an interaction between next 
nearest neighboring chains was included in the model. 

Recently, Ca3Co206 has drawn considerable attention 
to a step- like magnetization curve [f| H, 0, [1, H|. The 
number of the steps in the curve depends strongly on a 
sweep rate of the external magnetic field and tempera- 
ture [1, 0, 0- Two steps become apparent in the tem- 
perature range from 12 K to 24 K [9j. The first step 
takes place at the zero magnetic field. Then the mag- 
netic moment remains constant at about 1/3 of the full 
magnetization up to the magnetic field of 3.6 T where the 
second step occurs to the fully magnetized FM state. At 
least four equidistant steps are clearly visible below 12 
K at a very low sweep rate. They are accompanied by a 
sizeable hysteresis. Similar phenomena were observed in 
other spin-chain compounds, e.g. CaaCoRhOe [TcHTlj. 

The structure of CaaC^Og consists of C02O6 chains 
running along the c axis. The Ca ions are situated be- 
tween them. The chains are made up of alternating, face- 
sharing CoC>6 trigonal prisms and C0O6 octahedra. The 
crystalline electric field splits the energy level of Co 3+ 
ions into the high-spin (S — 2) and low-spin (S = 0) 
states. The Co 3+ ions situated in the trigonal environ- 
ment (Col) are in the high-spin state and the octahe- 
dral Co sites (Coll) occurs in the low-spin state. In the 
last case the energy difference between the low-spin and 



'Electronic address: kudasov@ntc.vniief.ru 



high-spin states is very small and a tiny fraction of Coll 
sites is reported to be in the high-spin state. The crys- 
talline electric field leads also to a very strong Ising-like 
anisotropy at the Col sites. The chains form triangu- 
lar lattice in the ab plane that is perpendicular to the 
chains. An in-chain exchange interaction between mag- 
netic Col ions through the octahedra with non-magnetic 
Coll ions is ferromagnetic (FM). The parameter of the 
FM in-chain coupling ( Jp) was found from the magnetic 
susceptibility at high temperatures [l2j , specific heat 0] , 
and theoretical calculations [l3[. These estimations are 
in a reasonable agreement with each other (Jp w 25K). 
The interchain interaction is antiferromagnetic (AFM) 
and much weaker than the in-chain one. It should be 
mentioned that the topology of the magnetic sublattice in 
CaaCo2 06 is much more complex than that in CSC0CI3, 
e.g. there exist helical paths [13| . A partial AFM order of 
chains appears at Tci—24 K. A weak feature concerned 
most probably with a transition in a new interchain order 
was also observed at around Tq2 ~ 12 K. This scenario 
of magnetic interactions in Ca3Co20g is consistent with 
results of x-ray photoemission spectroscopy [l4| , neutron 
scattering |15l|l6|. magnetization and specific heat mea- 
surements B S E3 j nuclear magnetic resonance 
(l8| . and theoretical calculations of indirect interactions 
between Col sites 

In the previous paper 19] I proposed a model to in- 
terpret the step-like magnetization curve of Ca3Co206 
at very low temperatures. It was based on two assump- 
tions, (i) The full in-chain FM order was supposed (rigid 
chains), i.e. only two state of a chain (spin-up and 
spin-down) were considered, (ii) Since Ca3Co206 and 
related compounds demonstrate very slow dynamics of 
chain re-orientation, the system was assumed to evolve 
over metastable states with a variation of the external 
magnetic field. A single-flip technique [20( was used to 
determine the metastable states as local minima. The 
ground state of 2D Ising model on the triangular lattice 
[2l| was taken as the initial state at the zero magnetic 
field. This approach allowed explaining the main features 
of the low-temperature step-like magnetization curve. At 
the same time, it is obvious that the first assumption is 
valid for the zero temperature only. In the present paper 
the rigid-chain model [l9j is extended to low non-zero 
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temperatures to clarify the origin of the transition from 
the low- to high-temperature magnetic structure at 12 K. 

Let us start with the rigid-chain model at the zero 
temperature. The ground state is described by the Ising 
model on the triangular lattice and can be represented 
through a sequence of approximations [ljj [2l[ . The hon- 
eycomb structure shown in the left panel of FigQ] is the 
first one. It can be improved by including an isolated tri- 
pod configuration (the middle panel of Fig. [1]), conjugate 
structures of two tripods, three tripods in a star config- 
uration (the right panel of Fig. [I}, and etc. In all the 
structures there exist chains with the zero effective field, 
that is, with three spin-up and three spin-down near- 
est neighbors. When the temperature is non-zero these 
chains should be in a disordered state. This fact breaks 
the first assumption of the rigid chain model. In a real 
crystal the chains are of a large but finite length. The 
upper limit of the number of atoms in a chain hardly ex- 
ceeds L = 10 3 because it is difficult to obtain the purity 
of a complex compound better than 99.9 %. In addition, 
taking into account structural defects we obtain L = 10 2 
as a reasonable estimation of the average chain length. 
Fortunately, it is shown below that the final results de- 
pend very weakly on the chain length. 

The partition function of the Ising chain at the zero 
effective field has a well-known form (e.g. [22T |) 



where 



2 L cosh^ 1 A 



(1) 



where K = Jp/T, T is the temperature, and Jf > 
for the FM in-chain coupling. A two-spin correlation 
function determined from Eq.([T]) is T(r) = tanh r K where 
r is the distance. Therefore, if T(L) « 1 there is a large 
probability of the full FM in-chain order. This condition 
can be written as 



K» iln(L). 



(2) 



It gives the crossover temperature at about 11 K that is 
close to the Tci- Above this temperature T(L) << 1 and 
we can't consider the chains at the zero effective field as 
rigid ones because they are in the disordered state. 

The most part of chains is under non-zero effective 
field. For i-th chain it can be written as 



<<T- > 



B 



(3) 



where j(i) denotes summation over the nearest-neighbors 
of the i-th chain, J a < is the parameter of the AFM 
interchain coupling, < cr| > is the average magnetic mo- 
ment of the j-th chain, and B is the external magnetic 
field. Under the effective field the partition function can 
be found exactly for the Ising ring (a chain with peri- 
odic boundary conditions) by means of a transfer matrix 
technique [22| 



Zc(H) 
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and Hi = h^/T. This equation gives the partition func- 
tion of the i-th single chain. For the sake of simplicity 
we have omitted the low indexes here. Due to the peri- 
odic boundary conditions Eq. (j4|) gives a different limit 
at H — ► as compared to that of Eq. ([I]) . However, both 
the expressions converge with increase of L. For very 
weak effective fields, i.e. sinh(iJ) ~ H « exp(— 2K), 
we obtain from Eq. (HJ the same condition for the FM 
in-chain order given by Eq. @. If the effective field is 
not extremely small, i.e. sinhi? >> exp(— 2K), the full 
FM order appears when tanh(iJ) ~ 2H » L^ 1 . We see 
that in any case a small effective field is sufficient for the 
crossover to the FM in-chain order. 

The foregoing discussion shows that at low tempera- 
tures we can assume the full FM in-chain order in all the 
chains with the exception of the chains with the zero or 
very small effective fields. Returning to FigfT] the state 
of the grey chains is rigid or disordered depending on the 
temperature. In fact this assumption is valid in a wide 
temperature range and broken only in a narrow region 
close to Tqi- Then the partition function of the honey- 
comb structure at the zero external magnetic field can be 
represented as 



Z hr = Z?/ 3 e- E °/ T 



(5) 



where N is the number of chains and Eq is the ground 
state energy of the honeycomb structure. At the zero 
temperature this expression gives the entropy She — 
(A/3) In 2 w 0.2317V. The tripod configuration appears 
at this limit with the probability 1/12 per chain and the 
entropy increases up to S t = (5JV/12) In 2 « 0.289iV 
At non-zero temperatures the tripod appears if the three 
chains in the centers of adjacent hexagons are in an or- 
dered state, the spin-down state (white circles) in the case 
shown in the middle panel of Fig. |T]). Thus, the prob- 
ability of tripod configuration per chain can be written 
as 



Pt = 



,K(L-l) 
Z r . 



= [2(1 
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Taking into account the the tripod configurations the 
partition function of the triangular lattice of the chains 
is given by 



z _ z N[l/3+P t ] e -E /T 



(7) 



It should be mentioned that the tripods do not exhaust 
all the configurations in the ground state of the 2D Ising 
model on the triangular lattice but they give a good ap- 
proximation as one can see comparing St with the exact 
value of the entropy of the triangular lattice S = 0.3231 A 

From Eq. one can easily calculate the free en- 
ergy and thermodynamic properties of the system. At 
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FIG. 1: The magnetic structures of the triangular lattice of chains: (left) the honeycomb structure, (middle) the isolated tripod 
configuration, and (right) the three tripod connected in the star arrangement. The black and white circles are spin-up and 
spin-down states, correspondingly. The gray circles can be either spin-up or spin-down states in the low-temperature phase or 
disordered chains in the high-temperature phase. 



the zero temperature this equation leads to the entropy 
(5iV/12)ln2. While the temperature grows Pt decreases 
and, when [1 + exp(— 2K)] L ~ 1 >> 1, the second term in 
the square brackets of Eq. ([7|) becomes negligible. Then 
Eq. ([7]) coincides with the partition function of the honey- 
comb structure ([5]). The crossover from the low- to high 
temperature regime occurs when exp(— 2K)(L - 1) « 1. 
This gives an estimation of the crossover temperature 
Tc2 ~H K very close to the experimental value. We 
also see that the dependence of the crossover tempera- 
ture on the chain length is logarithmic and, therefore, 
rather weak. 

Now we have a complete picture of the phase diagram 
in Ca3Co2C>6. At very low temperatures, i.e. much be- 
low Tc2, all the chains are rigid, that is, with the full 
FM in-chain order. The interchain magnetic order is 
the same as in the Ising model on the triangular lat- 
tice. While temperature increases the in-chain FM order 
collapses in the chains that are under the zero effective 
field. Whereas the entropy at the zero temperature is due 
to various configurations on the 2D triangular lattice, at 
higher temperature the entropy is accumulated mainly in 
disordered chains. That is why, the statistical weight of 
configurations with a larger number of disordered chains 
quickly grows with temperature. Since, as one can see 
from Fig. [TJ the honeycomb structure have the largest 
number of disordered chains, it dominates above Tq2 and 
the tripod configurations are suppressed. From this con- 
sideration we also see that the crossover at Tc2 is not a 
real phase transition. That is why, there is no peak or 
other features observed in the temperature dependence 
of the specific heat at the crossover temperature § . On 
the other hand, the crossover is rather sharp because the 
statistical weight of the tripod configurations depends on 
the temperature exponentially. 

Our interpretation of the magnetic structure transfor- 
mation agrees with the experimental observations. The 
two steps that are observed in the magnetization curve 
above T C2 1 together with the neutron scattering data 
[l6| indicate clearly the honeycomb magnetic structure. 
Using Eq. (|U) we calculated the magnetization curves at 
high temperature for the honeycomb structure as a initial 



•a i-o 



I 0.8' 
o 

¥ 0.6 
S 

I 0.4- 



tofj 
cs 



0.0-K 





i 


































A 








































' 1 


■ 





0.0 



1.2 2.4 
B(T) 



3.6 



FIG. 2: The magnetic moment as a function of the magnetic 
field. The solid line is the magnetization curve at a very 
low temperature [l9| . the dash line 2 and dash-dot line 3 are 
calculated for the honeycomb structure at T =12.5 K and 
T =20 K correspondingly. 



state (curves 2 and 3 in Fig. [J) . They fit the experimental 
curves [9(. This means that the smoothing of the steps 
at high temperatures is due to the magnetization of the 
disordered chains. On the contrary, the four steps cal- 
culated at the low temperature (Fig. [5]) are very sharp 
and some smoothing observed experimentally [3] stems 
most probably from the slow dynamics of the chains re- 
orientation. 

The magnetic phase diagram determined in thepresent 
work at H = is quite different from that of Ref . [3] . This 
discrepancy appears because of two reasons mainly, (i) 
As was mentioned above the structure of the magnetic 
sublattice in Ca3Co20fi is more complex as compared to 
CSC0CI3. (ii) In Ref. Q the in-chain exchange parame- 
ters was assumed to be equal to or even less then the in- 
terchain one. This gives rise to the randomly modulated 
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phase (RMP) in a wide temperature region. Fluctua- 
tions in the RMP develop along both a6-plane and c-axis. 
In Ca3Co206 the in-chain exchange parameter is much 
stronger than the interchain one. The RMP is suppressed 
in this case and the rigid honeycomb structure is stable 
up to high temperatures. Well above the crossover tem- 
perature our approach coincides with Shiba's mean-field 
model [4|. However, it can't be applied to the crossover 
region because all the phase boundaries in this model 
are phase transitions in contrast to the results obtained 
above and experimental data. 

The consequences of the results obtained are quite gen- 
eral. We can formulate universal rules for the phase dia- 
gram of frustrated Ising chain compounds, (i) The low- 
temperature magnetic phase is formed by rigid chains 
and its ground state is determined by the ground state 
of the 2D Ising model. It should be mentioned that in 



case of Mekata's theory of CsCoCbj the ground state 
of the 2D Ising model is different from that considered 
in the present paper because of the additional interac- 
tion between next-nearest- neighbor chains, (ii) In the 
high-temperature magnetic structure the number of dis- 
ordered chains at the zero effective field that should be 
maximal to increase the entropy and lower the free en- 
ergy. That is why, at high temperatures the magnetic 
structure of the triangular lattice of Ising chains tends 
to the honeycomb one. These rules are valid for the 
FM and AFM Ising chain systems Ca 3 Co 2 6 , CsCoCl 3 , 
CsCoBr3 as well Ising chain compounds with distortions, 
e.g. TlCoC13 [13 and RbCoBr3 H. 

I am grateful to V. Uzdin for fruitful discussions. This 
study was partially supported by INTAS grant (03-51- 
4778) "Hierarchy of scales in magnetic nanostructures" . 
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